Generation of non-classical photon states in superconducting quantum metamaterials 
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We report a theoretical study of diverse non-classical photon states that can be realized in 
superconducting quantum metamaterials. As a particular example of superconducting quantum 
metamaterials an array of SQUIDs incorporated in a low-dissipative transmission line (resonant 
cavity) will be studied. This system will be modeled as a set of two-levels systems (qubits) strongly 
interacting with resonant cavity photons. We predict and analyze a second-order phase transition 
between an incoherent (the high-temperature phase) and coherent (the low-temperatures phase) 
states of photons. In equilibrium state the partition function Z of the electromagnetic field (EF) in 
the cavity is determined by the effective action Sef f{P{T)} that, in turn, depends on imaginary-time 
dependent momentum of photon field P{t). We show that the order parameter of this phase 
transition is the Po{t) minimizing the effective action of a whole system. In the incoherent 
state the order parameter Poir) ~ but at low temperatures we obtain various coherent states 
characterized by non-zero values of Po{t)- This phase transition in many aspects resembles the 
Peierls metal-insulator and the metal-superconductor phase transitions. The critical temperature of 
such phase transition T* is determined by the energy splitting of two- level systems A, a number of 
SQUIDs in the array , and the strength of the interaction rj between SQUIDs and photons in cavity. 

(submitted to Superconductor Science and Technology) 

PACS numbers: 42.50.-p,74.81.Fa,74.50.-|-r 
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I. INTRODUCTION 

A great attention is devoted to a theoretical and ex- 
perimental study of novel quantum metamaterials^"— . 
These systems consists of a large amount of solid-states 
elements (qubits), i.e. two-levels systems showing di- 
verse coherent quantum phenomena, e.g. quantum beat- 
ing, microwave induced Rabi oscillations, Ramsey fringes 
etc4. To obtain such coherent quantum-mechanical be- 
haviour in single qubits, the dissipation and decoherence 
have to be small enougbi. Moreover, in order to observe 
novel collective coherent quantum effects in metamate- 
rials a strong long-range interaction between single ele- 
ments has to be provided by surrounding media. 

Various superconducting systems, e.g. arrays of 
Josephson junctions, RF SQUIDs, many-junctions super- 
conducting quantum interferometers, just to name a few, 
incorporating in low-dissipative (superconducting) trans- 
mission line are extremely suitable in order to realize 
such quantum metamaterials. Indeed, dc-biased Joseph- 
son junctions^- or diverse SQUIDsii^ subject to an ex- 
ternally applied magnetic field have been established re- 
cently as qubits in which various coherent quantum ef- 
fects have been observed. A strong long-range interaction 
between well-separated qubits is provided by a transmis- 
sion line through emission (absorption) of virtual pho- 
tons. This type of interaction was proposed in Refsi^"— 
and realized in experiments with single qubits incorpo- 
rated in a resonatorii^ii^. It has been shown that a strong 
interaction between well-separated qubits results in an 
enhancement of quamtum-mechanical tunnelingi^ii^iii 
and suppression of decoherence induced by a spread of 



parameters of qubitsiS. The measurements of frequency 
dependent transmission (reflection) coefficient of electro- 
magnetic field (EF) propagating through the transmis- 
sion line provides a convenient method to observe coher- 
ent quantum phenomena in such metamaterial o^^'^" . 

On other hand a strong interaction of qubits with EF 
can result in different states of photons in the cavity. In- 
deed, in the absence of interaction with qubits the unique 
photon state of the cavity is an incoherent, chaotic state 
of photons characterized by a well-known Planck distri- 
bution, i.e. < E >^ 0, < {Ef cx [e^'^ - l]-^ >, 
where E is the electric field of radiation, < ... > is the 
quantum-mechanical average, and wq is the frequency of 
the cavity mode. A strong interaction of EF with qubits 
leads to the effective enhancement of energy levels dif- 
ference of qubits that, in turn, changes EF in the cavity. 
Thus, one can expect that in the resonant cavity strongly 
interacting with an array of qubits the different states of 
photons can be observed. In this paper we analyze pos- 
sible quantum-mechanical states of photons emerging in 
the resonant cavity strongly interacting with an array of 
qubits. Notice here, that a similar analysis of the photon 
states of the cavity interacting with an unbiased array 
of Josephson junctions has been done in-- in order to 
explain a strong radiation from a 2d-array of Josephson 
junctions observed in Ref.— . 

The paper is organized as follows: in Section II we 
present a model of an array of SQUIDs incorporating in a 
low-dissipative transmission line, and elaborate the clas- 
sical description of this system, i.e. the dynamic equa- 
tions of motion and the Lagrange function, in Section 
III a complete quantum-mechanical description of this 
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FIG. 1: The schematic of an array of RF SQUIDs incorpo- 
rated in a transmission line. 



model is provided. In section IV using a great similar- 
ity with well-known phase transitions, e.g. the metal- 
ferromagnelj^^, metal-superconductor— and the Peierls 
metal-insulator— transitions, we predict and analyze a 
second-order phase transition between the incoherent, 
chaotic state (the high-temperature phase) of photons 
and diverse coherent non-classical photon states (the low- 
temperature phase). The Section V provides discussion 
and conclusion. 



II. MODEL AND CLASSICAL DESCRIPTION 
OF SUPERCONDUCTING METAMATERIALS 

As a particular example we consider here a system 
of RF SQUIDs incorporated in a low-dissipative trans- 
mission line. An each RF SQUID is characterized by a 
Josephson phase ipi = 27r$i/(/)o, where $i is the total flux 
in the superconducting loop of a SQUID, and (fio is the 
flux quantum. An application of dc-magnetic field char- 
acterized by ^ext allows one to tune the potential relief of 
a Josephson phase cpi from a single well up to a double- 
well potential. The set of RF SQUIDs is incorporated in 
a linear transmission line. The transmission line is char- 
acterized by two parameters Lq and Co, the inductance 
and capacitance per unit length, accordingly. We also in- 
troduce the voltage V{x) and current I{x) distributions, 
where x is the coordinate along a transmission line. The 
inductive coupling, M — tjLq, provides an interaction be- 
tween RF SQUIDs and transmission line. The schematic 
of a system is presented in Fig. 1. 



A. Classical Equations of Motion: linear 
transmission line 



We start with the classical dynamic equations for a 
linear transmission line. It is 



dV{x,t) dl{x,t) 



dx 



dt 



and 



dl{x,t) ^ dV{x,t) 

Cq- 



dx 



dt 



(1) 



(2) 



These two equations are rewritten as 
d^Iix,t) 1 d^I{x,t) 



dx"^ 



9t2 



(3) 



where LqCq = l/c^. The electromagnetic standing waves 
can occur in this 1-d cavity resonator. The wave vec- 
tors are determined by standard boundary conditions: 
kn — 7rn/L, where L is the size of a transmission line, 

n = 1,2 We will consider a transmission line with an 

extremely high quality factor, which was routinely ob- 
tained in superconducting transmission lines, and there- 
fore, the only one wave vector will be important in the 
dynamics of coupled RF SQUIDs and EWs of the cavity. 



B. Classical Equations of Motion: an individual 
SQUID 

The classical dynamic equations for an individual 
SQUID are written as 



IT ■ / N oi dip 1 (fip 
Irf/Ic = sm((^) + —— + — — 



(4) 



where 1^ is the critical current of a Josephson junc- 
tion, \/a^ is the McCumber parameter characterizing 
the dissipation of SQUID, is the plasma frequency of a 
Josephson junction. On other hand the Josephson phase 
in a SQUID loop is satisfied to the following equation as 



V = fext - PlIrf/Ic , 



(5) 



where /3l is the inductive (dimensionless) parameter of 
the SQUID, fext corresponds to the sum of the externally 
applied dc-magnetic field and ac-magnetic field induced 
by a current flowing along the transmission line. Thus, 
(ysf^t allows one to tune the potential relief of the Joseph- 
son phase, and (p'^^^ cx / provides a coupling between 
the RF SQUID and a transmission line. 



C. Classical Equation of Motion: coupled 
transmission line and SQUIDS 

Inductive coupling between RF SQUIDs and transmis- 
sion line results into a particular change of classical equa- 
tions of motion: first, the Eq. ([3]) changes to 



d^I{x,t) 1 d^Iix,t) 



dx"^ 



dt^ 



E 



c2 at2 



(6) 



where rj is the parameter characterizing a mutual induc- 
tance of RF SQUIDs and the transmission line; secondly, 
=(ryLo/</'o)/(a;.,i). 
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D. Lagrangian of a Superconducting Metamaterial and 



The classical equations of motion can be derived also 
from the Lagrangian of a whole system 
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2Co ' dx 



- "Pelt + ivLo/MQi^,t) 



(1 - COS<^i)7 



(7) 



where Q{x,t) is the charge variable characterizing a 
transmission line, and Ej — hlc/{2e) is the energy of 
Josephson junction. To simplify this expression we con- 
sider the standing EWs with a single wave vector fc„ 
only. In this case the charge variable Q{x,t) has a form: 
Q(x,t) = Q (t) cos{knx) . Substituting this expression in 
Eq. d?]) we obtain 



where 

Lph = rn 
and 



L — Lph + LsQuiD + Li 



= Lo/2, (8) 



LSQUID = EJJ2:^[^^f-^[V^-^^xt?-{^-C0S^^), 

(9) 



and 



Lint = - '^{r]Lo/^o)-T^y^i cos{knXi)Q{t), (10) 
i 

where the coordinates of RF SQUIDs along a trans- 

mission line. 



P{t) 



Using these new variables the photon Hamiltonian is 
written as 



B. RF SQUID Hamiltonian 



(12) 



We consider the macroscopic quantum dynamics of RF 
SQUID when the potential energy has a form of double 
well potential. In this case the Hamiltonian of a system 
of isolated RF SQUIDs is written as 



H SQUID = ^ -[AjCTa; + ejCT^] , 



(13) 



where Uz and Ux are standard Pauli matrices, Aj and 
ti are the matrix element of tunneling and the energy 
difference between two potential wells, accordingly. No- 
tice here, that such a Hamiltonian can be used also for 
more complex qubits, e.g. phase qubits, flux qubits etc, 
where parameters and are determined by the phys- 
ical properties of corresponding qubits. For qubits based 
on RF SQUIDs we obtain 



A cx fiujp{l 



1//3l) 



and 



cx 7rV6(l-l//3L)'/'(</'ext/</'o-l). 



(15) 



Moreover, the parameters A^ and Ci can fluctuate from 
one qubit to other one. 



Interaction Hamiltonian 



III. QUANTUM DESCRIPTION OF 
SUPERCONDUCTING QUANTUM 
METAMATERIAL S 

A. Photon Hamiltonian 

Introducing the "momentum " of a photon field P{t) as 
P = {Lo/2)Q{t) we obtain the Hamiltonian of a photon 
field in the following form: 



^ph — 



(11) 



La ' 4 

Next we introduce the operators of boson field b and b 



The equilibrium dynamics of a Josephson phase in 
imaginary time representation can be presented as rare 
jumps (the instanton type of solution) between two equi- 
librium positions^. Using this property we obtain the 
interaction Hamiltonian as follows: 



H^nt = « > ^i(Jz (6 - 6+) , C» = 7 (16) 

~^ PL 00 

where Sip is the phase difference between two equilibrium 
positions of a Josephson phase. 



D. Effective action 



Q{t) 



cknLo 



(6 + 6+) 



In order to study the various photon states arising in 
superconducting quantum metamaterials we write the 
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partition function Z in the form of functional integral 
as 

Z = J DQD{ip,}expi~S{Q,ip,}), (17) 

where S is the action of EF interacting with an array of 
two- levels systems. Integrating the Eq. ([T7|) over {(pi} we 
obtain the effective action Sef / in the following form~: 



SefAQir)] 



fi Jo 



h/(kBT) 



dr- 



/jsT^ln 



cosh - 



(18) 



where Q!i{(5(T)} are the positive Floquet eigenvalues of 
arrays of two-levels systems in the presence of periodic in 
imaginary time potential Q{t). Notice here, that in the 
absence of interaction with the field Q{t), i.e. as = 0, 
the Floquet eigenvalues are ai(0) = l^f + ef, and the 
minimum of Seff[Q{T)] occurs as Q = 0. 



IV. PHASE TRANSITIONS IN STATES OF 
PHOTONS 

Since the interaction term in the Lagrange function 
Lint depends on r-derivative of Q, i.e. Q, we introduce a 
new variable of the momentum of a photon field P{t) = 
mQ and rewrite the effective action Sef / as 



1 rh/ikBT) , 



UJn 



fcsT^ln 



cosh 



(19) 



Moreover, the Floquet eigenvalues ai{P) are determined 
from the following equation; 



{dr + H'p)i>, = 



and the Hamiltonian Hp is 



El + mP 

A,, 



A. 

-(e. + f/zP) 



(20) 



(21) 



where the parameter determining the interaction be- 
tween cavity modes and two-levels systems ry,; = 
''^*ofe cos(fc„a;i). 

Next, we obtain the periodic in imaginary time repre- 
sentation function Pq{t) minimizing the effective action, 
as a solution of the following equation: 



Pa ^ 1 p . L 

m mun '—^ oP 



n{Po} 



(22) 



A. Classical Second Order Phase Transition 

First we consider a particular case as the momen- 
tum of a photon field Pq does not depend on r. 
The Floquet eigenvalues are determined as ai{P) — 
\J (ci + ?Ti-Po)^ + Af and the equation for Pq reads as 



_ - (fj + T?iPo) 



tanh 



(23) 

In a simplest case as = and fji = 77 we obtain self- 
consistent equation 



Po 



1 



E 



1 



._tanh ^ ' " = 



(24) 

At high temperatures, T > T*, this equation has a single 
solution Pq — 0. Such a high-temperature phase corre- 
sponds to the incoherent, chaotic state of a photon field. 
However, at low temperatures, T < T*, the Eq. (|24|) has 
two non-zero solutions of ±Poi and t/ie coherent state 
of a photon field can be realized. Therefore, we obtain 
a second-order type of phase transition in the states of 
photon field interacting with a set of two-levels systems. 
The critical temperature T* depends on the distribution 
of energy differences of two- levels system A^ . If such dis- 
tribution is a narrow one, i.e. A^ ~ A, we obtain the 
value of T* as 



mfj'^N 



(25) 



where is a total number of two-levels systems. Such 
a phase transition occurs only if Aq < ksT*. At low 
temperatures Pq reaches the maximum value of Pq = 
-iimfjN . As one can see the maximum value of Pq is 
proportional to N , and this dependence also indicates the 
presence of the coherent state. The dependence of Pq{T) 
for this case is shown in Fig. 2 (red solid line). Notice here, 
that this case resembles a well-known metal-ferromagnet 
phase transition^,. 

In the opposite case of a wide distribution of Ai, e.g. 
from zero to upper cut-off Aq, the critical temperature 
is determined as 



(26) 



This phase transition occurs if Aq > mrfN. At low 
temperatures the Pq reaches the maximum value of 
I Pol ~ ksT^/fj. This phase transition resembles 
the superconductor- normal metal^l and/or Peierls metal- 
insulator^^ transitions. 



B. Quantum Phase Transitions 

Here we consider a phase transition into a peculiar 
"quantum ordered" state representing a quantum inter- 
ference between the two semi-classical states ±Po of the 
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FIG. 2: The classical and quantum phase transitions in the 
state of photons. The temperature dependencies of the mo- 
mentum of photon field Po{T) are shown: r-independent Po 
(red line) and Po(t) oc Pq (T) sin(Po (T)r) (blue line). The 
case of a narrow distribution of A,; ~ Aq is shown. The 
parameters = 2QK, and Aq = 'iK were used. 



photon field, which are inversion symmetry related solu- 
tions of the self-consistency equation Eq. ([M| . Each of 
the ±Po states separately describes a particular coher- 
ent state of the photon field in our model. The quantum 
ordered state was predicted in22 for a system of electron- 
hole pairs coupled to a semi-classical spin-density wave 
fluctuations. In our present model it is described by 
the amplitude of the semi-classical photon field Po{t)^ 
or the " quantum order parameter" , which is the periodic 
in imaginary time solution of Eq. (|22p . Again we con- 
sider a simplest case as = and f/i = t) and apply 
analytic solution found in^^: 

fjPoiT) = AnKTkisn {inKTr; h) , K ^ K{ki) (27) 

where sn(r, fci) is the Jacobi snoidal elliptic function, pe- 
riodic in T with period l/(nT), n = 1,2, K{ki) is el- 
liptic integral of the first kind, positive integer n and real 
parameter < fci < 1 are found by minimizing the Eu- 
clidian action S^f / given in Eq. (|19p . Here we merely de- 
scribe a single-harmonic limit fci — > of solution Eq. (P?]) . 
In this limit expression for i-b(''') in Eq. (P7| turns into: 



f]Po{T) w 2'KnTkisin {2'KnTT) . 



(28) 



Simultaneously, it was shown in^^, that solution in the 
form of Ea. (l?7l) leads to the following spectrum of the 
Floquet eigenvalues ai{P) found from Ea.pO| : 



ai 



2A, 



1 ~ fc^ + A^ 
1 + A? 



1/2 



1 + A' 



(29) 



The latter relation between parameters fci , fc is known as 
Landen transformation^-. The Jacobi's function M(r) = 
M(t + 1/nT) from Eq. ([27l) turns the generic self- 
cosistency equation (|22p into algebraic equation for pa- 
rameters k, v?'^ : 



E 



tanh ■ 



A, 



1 



{(Af + l)(A| + fc'2)}i/2 



(31) 



It is not hard to see, that in the limit fc' — >■ Eq. (PT|) 
transforms into classical mean-field self-consistency Eq. 
p4| . While in the limit fc' — >■ 1 the self-consistency equa- 
tion Eq. as it follows from Eq.dUl) and Eq. ljHO)) 
turns into the following equation: 



E 



tanh ■ 



A,, 



knT 



1 



Af + {^nPy 



mq 



n2 



(32) 



Now, after a comparison with the Eq. ([Ml) for the classical 
photon condensate Pq it is possible to conclude, that in 
the case of narrow distribution of energy differences of the 
two-levels systems A^ w Aq and strong coupling constant 
to the electromagnetic field: Aq < raff^N , the quantum 
ordered phase (QOP) of the photon field occurs below 
the temperature: 



rpQOP 



1/3 



1 

ks 



(33) 



Since the amplitude of PaiT) is proportional to AnKTki 
in accord with Eq. (P7)) . this result suggests that the 
number of photons condensed into quantum ordered 
phase is cx N'^/'^^ where TV is the number of the two-lewel 
systems. Hence, this dependence indicates the presence 
of the coherent state also in the quantum ordered phase, 
but with less strong entangling than in the classical pho- 
ton condensate described by Eq. p4)). In the opposite 
case of a wide distribution of A^ and weak coupling con- 
stant Aq > raff^N we find transition temperature similar 
to the classical mean-filed case Eq. ([25]) : 



(34) 



In both cases, when temperature lowers well below T^^^ 
the increase of the integer number n cx 1 /T of the oscil- 
lations of the quantum order parameter P(t) along the 
imaginary time-axis t within the Euclidian space tempo- 
ral interval [0,1/r] keeps the QOP amplitude 4nKTki 
finite and non- vanishing up to T = OK state acording to 
Eq. (j27p . Thus, we can conclude that the quantum phase 
transition in the photon states is the first-order type of 
phase transition. 



Here H(m, k) is elliptic integral of the third kind, and 
besides: 



A,; 



2TnK{ky 



, k = 2Vfci/(l + fci); fc'^ = 1 - fc" (30) 



V. DISCUSSION AND CONCLUSIONS 

In previous Section we obtained classical and quantum 
second-order phase transitions in the states of photons 
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that can be emerged in a resonant cavity strongly in- 
teracting with an array of two-levels systems. For the 
classical phase transition we obtain that at low tempera- 
tures, T < T*, the incoherent photon state (with Pq = 0) 
becomes unstable. Two novel photon states character- 
ized by non-zero values of the momentum of photon field 
±Pq(T) can appear at low temperatures. These photon 
states are the coherent photon states well-known in the 
quantum optics^S. The average number of photons in 
these states n = Pg /(2?7ifia;o). The temperature depen- 
dence of Po{T) is determined by Eq. ([Ml) , and it is shown 
in the Fig. 2 (red solid line). The critical temperature 
T* depends strongly on the distribution of energy dif- 
ferences Ai, a number of qubits TV, and the strength of 
interaction f/. Thus, for reasonable values of parameters 
T* varies from lOOmK up to 50K. The value of T* is 
determined by Eqs. ((25|) and (|26|) . This phase transition 
is similar to a well-known metal-fcrromagnet transition, 
and the momentum of photon field ±Po {T) corresponds 
to the " Weiss effective magnetic field" in a theory of fer- 
romagnetic materials^. 

Notice here, that these two coherent photon states cor- 
responding to the -l-Po(T) and — Po(r) values are the 
degenerate ones. In this case the quantum beating be- 
tween these two states can be observed in the system. 
The quantum beating between two photon states results 
in a splitting of resonant frequencies of the cavity, i.e. 
ujres — ^0 ^ Aw. The splitting Auj is obtained as fol- 
lows: the two stable photon states are separated by the 
effective potential barrier AU ~ mN'^fj^ (for T — Q), 
and therefore. Aw ~ wq exp[— AJ7/(?ia;o)]. The effect of 
the splitting of resonant frequencies is a consequence of 
the degeneracy of the photon states. Such a degeneracy 
can be lifted by application of an external magnetic field 
allowing one to realize non-symmetric double- well poten- 
tial for RF SQUIDs, i.e. as e,; ^ 0. In this case a single 
coherent state of photons with Pq cx ti emerges in 
the cavity. 

We obtained also that different metastable states of 



photons can be obtained in this system. In these states 
there is no net classical photon condensate, but there 
exists macroscopic quantum condensate (with amplitude 
of the photon momentum P c>c iV^/^, where N is the 
number of two-level systems), that has zero mean value 
of the electromagnetic field. These states appear as a 
result of a first-order phase transition. 

In conclusion we have shown that superconducting 
quantum metamaterials can support the diverse non- 
classical photon states. As a particular example of such 
metamaterial we considered an array of RF SQUIDs 
incorporated in a low-dissipative resonant cavity. We 
mapped this system to a set of two- levels systems (qubits) 
strongly interacting with photons of the cavity. By mak- 
ing use of a complete quantum-mechanical description 
of such a system we obtained that at high temperatures, 
T > T* the incoherent chaotic state of photons is a stable 
one. At low temperatures T < T* a. large amount of dif- 
ferent photon states emerges in the cavity. These photon 
states appear as a result of specific classical (the second- 
order type) or quantum (the first-order type) of phase 
transition. The physical origin of such phase transitions 
is a following: a strong interaction of EF with two-levels 
systems leads to the effective enhancement of energy lev- 
els difference of qubits that, in turn, changes the EF in 
the cavity. The order parameter of phase transitions is 
the T-dependent momentum of photon field Po(''"). In the 
case of the classical phase transition as Pq (t) — const the 
coherent photon states and the quantum superposition 
of two coherent photon states can be obtained in the cav- 
ity. In the case of quantum phase transitions the differ- 
ent metastable photon states characterized by complex 
dependence of P(t) (see Eq. ([17])) also can be realized. 
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^ J. Plumridge, E. Clarke, R. Murray, and Ch. Phillips, Sol. 
State Comm,146 406 (2008). 

^ A. L. Rakhmanov, A. M. Zagoskin, S. Savelev, and F. Nori, 
Physical Review B 77 144507 (2008). 

^ J. Q. Quach, C-H. Su, A. M. Martin, A. D. Greentree, and 
L. C. L. HoUenberg, Optics Express 19 11018 (2012). 

* A. M. Zagoskin, Quantum Engineering: Theory and De- 
sign of Quantum Coherent Structures., Cambridge Univer- 
sity Press, Cambridge, 272311 (2011). 

^ J. Martinis, John, S. Nam, J. Aumentado, and C. Urbina, 
Phys. Rev. Lett. 89, 117901 (2002). 

^ J. Lisenfeld, A. Lukashenko, M. Ansmann, J. M. Martinis, 
and A. V. Ustinov Phys. Rev. Lett. 99, 170504 (2008). 
I. Chiorescu, Y. Nakamura, C. J. P. M. Harmans, and J. 
E. Mooij, Science 299, 1869 (2003). 

8 R. H. Koch, J. R. Rozen, G. A. Keefe, F. M. MiUiken, C. 
C. Tsuei, J. R. Kirtley, and D. P. DiVincenzo, Phys. Rev. 



Lett. 96, 127001 (2006). 
^ Al. Blais, R.-Sh. Huang, A. Wallraff, S. M. Girvin, and R. 
J. Schoelkopf Phys. Rev. A 69, 062320 (2004). 
S. H. W. van der Ploeg, A. Izmalkov, Alec Maassen van 
den Brink, U. Hiibner, M. Grajcar, E. Ilichev, H.-G. Meyer, 
and A. M. Zagoskin, Phys. Rev. Lett. 98, 057004 ( 2007). 
" M. V. Fistul and A. V. Ustinov, Phys. Rev. B 68, 132509 
(2003). 

^2 M. V. Fistul and A. V. Ustinov, Phys. Rev. B 75, 214506 
(2007). 

M. V. Fistul, Phys. Rev. B 75, 014502 (2007). 
" A. WaUraff, D. I. Schuster, A. Blais, L. Frunzio, R.-S. 
Huang, J. Majer, S. Kumar, S. M. Girvin, et al. Nature 
431, 162 (2004). 

I. Chiorescu , P. Bertet, K. Semba, Y. Nakamura, C. J. 

Harmans, and J. E. Mooij, Nature 431 159 ( 2004). 

X. Y. Jin, J. Lisenfeld, Y. Koval, A. Lukashenko, A. V. 



7 



Ustinov, and P. Miillcr Phys. Rev. Lett. 96, 177003 (2006). 
S. Savelev, A. L. Rakhmanov, and F. Nori, Phys. Rev. 
Lett. 98, 077002 (2007). 

S. Ghosh, W.H. Wang, F. M. Mendoza, R. C. Myers, X. Li, 
N. Samarth, A. C. Gossard, and D. D. Awschalom, Nature 
Materials 5, 267 (2006). 

A. A. AbdumaUkov, Jr., O. Astafiev, A. M. Zagoskin, Yu. 
A. Pashkin, Y. Nakamura, and J. S. Tsai Phys. Rev. Lett. 

104, 193601 (2010) 

M. Jerger, S. Poletto, P. Macha, U. Huebner, A. 
Lukashenko, E. Il'ichev, and A. V. Ustinov, Europhys. 
Lett. 96, 40012 (2011). 

J. Kent Harbaugh and D. Stroud, Phys. Rev. B 61, 14765 
(2000). 

P. Barbara, A. B. Cawthorne, S. V. Shitov, and C. J. Lobb 



Phys. Rev. Lett. 82, 1963 (1999). 

Ch. Kittcl, Introduction to Solid State Physics, 8th. ed, 
Wiley, N. Y. (2005) 

K. Machida and H. Nakanishi, Phys. Rev. B 30, 122 

(1984) . 

K. Machida and M. Fujita, Phys. Rev. B 30, 5284 (1984). 
S. Chakravarty and S. Kivelson, Phys. Rev. B 32, 76 

(1985) . 

S. I. Mukhin, J. Supercond. Nov. Magn 24, 1165 (2011). 
S. I. Mukhin, J. Supercond. Nov. Magn 22, 75 (2009). 
E. T. Watson and G. N. Witteker, Modern analysis, Cam- 
bridge University Press., Cambridge, 1952. 
R. Loudon, The Quantum Theory of Light, Oxford Univer- 
sity Press, N. Y. (1973). 



